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AN APPROXIMATE SOLUTION FOR THE FREE VIBRATIONS OF 
ROTATING UNIFORM CANTILEVER BEAMS 


David A. Peters 
Ames Research Center 
and 

U.S. Army Air Mobility Research and Development Laboratory 
Moffett Field, California 94035 


SUM4ARY 

Approximate solutions are obtained for the uncoupled frequencies and 
modes of rotating uniform cantilever beams. The frequency approximations 
for flap bending, lead>lag bending, and torsion are simple expressions having 
errors of less than a few percent over the entire frequency range. These 
expressions provide a simple way of determining the relations between mass 
and stiffness parameters and the resultant frequencies and mode shapes of 
rotating uniform beams. 


INTRODUCTION 

Many theoretical analyses of hingeless rotors are performed for the ideal 
case of uniform untwisted rotor blades. These analyses need to be performed 
over a wide range of blade stiffness parameters requiring calculation of a 
broad spectrum of rotating blade natural modes and frequencies. In such 
theoretical investigations, the investigator often needs to translate some 
desired set of natural frequencies into blade properties and natural mode 
shapes which correspond to those frequencies. Then the results can be com- 
pared with other theories as well as experiments. There is, therefore, a 
practical requirement for a simple and accurate means of relating the fre- 
quencies and modes to the mass and stiffness characteristics of rotating uni- 
form beams. 

The usual methods for calculating frequencies and modes of rotating beams, 
however, are not well suited to the particular computations required in this 
type of analysis. First, traditional energy and finite element methods yield 
poor convergence in the practical range of flapping frequenwies (and in the 
lower extremes of lead-lag and torsion frequencies). Thus, a moderate number 
of assumed modes (or finite elements) must be used. As pointed out in Reference 
1, this poor convergence necessitates an involved numerical solution for the 
frequencies and modes, which must be repeated for each separate blade configura- 
tion. Second, because the solutions are numerical, a trial and error process 
is required to determine the blade properties which will result in a desired 
set of natural frequencies. 

The purpose of this report is to provide simplified expressions for the 
frequencies and modes of rotating uniform beams. It will be shown that these 
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expresisons are accurate through the entire blade stiffness range yet elementary 
enough to overcome the disadvantages of purely numerical solutions. 

Equations of Motion 

The uncoupled equations ("uncoupled" implies: zero twist, coincident mass 

center and elastic axis, beam parallel to plane of rotation, and warp included 
only in calculation of J) for the free vibrations of rotating uniform beams 
are< 

nu - y ((1 - r^)u'l' - w^u « 0 (flapping) (la) 

nv**" “ 2 ~ - (u^ + l)v « 0 (lead-lag) (lb) 

-y 6" - j [(1 - r^)e']’ - (o)2 - l)e = O (torsion) (Ic) 

where u and v are the flapping and lead-lag displacements nondimensionalized 
on the blade radius R, r is the distance along the beam divided by R, u is 
the natural frequency of the vibration divided by the rotational frequency 
n, and n and y are stiffness parameters. 


n ® 


El 


y a GJ 

mJ22R2lc2 


The first term in each line of Eq. (1) represents the elastic stiffness of the 
blade in bending or torsion. The term in each equation represents the 
inertial force due to the harmonic motion. These terms dominate for large 
values of n and y. The remaining terms are due to the centrifugal tension, 
and they dominate (with u^) for small values of n and y. The boundary condi- 
tions corresponding to a cantilever beam are 

u(0) » u'(0) * u"(l) « u"*(l) * 0 (2a) 

v(0) « v'(0) » v"(l) = v'"(l) * 0 (2b) 

0(0) = e»(l) « 0 (2c) 

The r = 0 conditions are geometric boundary conditions which specify that the 
displacement (and slope in bending) must vanish at the root. The r » 1 con- 
ditions are natural boundary conditions which specify that the shears and 
moments vanish at the tip. Eq. (1) and Eq. (2) together form the uncoupled 
rotating beam equations. 
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Solution by Orthogonal Expansion 


Eq. (1) is often solved using an expansion in orthogonal functions which 
satisfy the geometric boundary conditions. With application of an energy 
method (e.g. Galerkin or Rayleigh-Ritz) , such expansions transform the charac- 
teristic differential equation into an algebraic eigenvalue determinant.^ 

When the known, exact nonrotating mode shapes 41 ,^ are chosen as the orthogonal 
functions, this determinant takes on the particularly simple form 




where are the nonrotating frequencies and 


0 



“ - r2)dr,A « 


(flapping) 
♦ 1 (lead- lag) 
- 1 (torsion) 


(3) 


(4) 


In practice, [D] is trancated at a specified number of rows and columns. Its 
eigenvalues and eigenvectors aj^j then give approximations for the rotating 
frequencies u^, Eq. (4), and for the rotating mode shapes. 
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n 

In Table I, formulas are given for the nonrotating modes 4 )., nonrotating 
frequencies and the Galerkin integrals A^j. Ihe formulas for nonrotating 
frequencies are specifically 

u^Nr ® (flap, lead- lag) 

2 ir^(2n - 1)^ ^ 

“nr ' — ^""4 ^ ^ (torsion) 

Whenever convenient, therefore, the stiffness parameters n and y may be 
expressed solely in terms of nonrotating frequencies and universal constants. 

Numerical Limitations 

In Fig. 1, the value of the lowest frequency of the eigenvalue determi- 
nant is presented for various modal truncations of the flapping equation.^ 

As decreases, the number of orthogonal functions N must be increased in 
order to retain a specified accuracy. In fact, when = 0, no value of N 
will result in the known exact value of the first frequency, u « 1.0. This 
failure at Uwj. - 0 can be explained in mathematical terms. When n ° 0, 






Eq. (la) is reduced from a fourth order equation to a singular second order 
equation, relaxing all but one of the boundary conditions of Eq. (2a). Thus, 
the orthogonal functions chosen to satisfy the boundary conditions of the 
fourth order equation are not applicable to the second order case. From a 
physical standpoint, the n " 0 modes (which have nonzero slopes at r » o) 
cannot be formed from a linear combination of nonrotating cantilever beam 
modes (which have identically zero slope at the root). 

Although this failure of the eigenvalue determinant at ■ 0 does not 
strictly hold in the range of practical flap frequencies (0.05 < < 0.50). 

such values of are close enough to zero to create convergence difficulties 
in the eigenvalue determinant. Fig. 1 reveals that ten or more orthogonal func 
tions may be required for satisfactory results. Hius, involved numerical cal- 
culations are required in order to find a.. Furthermore, the number of orthog- 
onal functions required in this procedure is also an unknown which must be 
determined by iteration. These same drawbacks also apply to the lead-lag and 
torsion equations (although to a much lesser degree) . It is the purpose of 
this paper to provide simple expressions for u, valid at all values of so 
that the above drawbacks can be overcome. The flapping equation will be con- 
sidered first, with the lead-lag solution following directly by replacing 
“flap ^ * “^LEAD-LAG torsion equation will be 

considered. 


SMALL STIFFNESS EXPANSION FOR BENDING 

The first step toward developing simple frequency expressions is the 
determination of the asymptotic behavior of u^(flapping) as (i.e. n) 
approaches zero. A similar phenomenon, the asymptotic behavior of a stat- 
ically loaded beam under tension as El approaches zero, is a classic problem 
in the theory of matched asymptotic expansions.*^ That solution indicates 
that as El approaches zero the beam takes on the shape of a loaded string 
(El • 0. nonzero slope near the supports) except in a small region adjacent 
to the supports where the boundary condition (zero slope) is fulfilled. 

Since the rotating beam flapping equation is similar in character to this 
classic example, a matched asymptotic expansion promises to provide the 
desired small n solution. 


Outer Expansion 

General solution . A dimensional analysis of Eq. (1) indicates that the 
basic length dimension (divided by R) of the flapping equation is Thus, 

an expansion in powers of is a logical starting point for the analysis. 
Expanding u and in 

u = u^ ♦ + nu^ + n^^^Uj + . . . (5a) 

0)2 a + . . . (5b) 

0 12 3 ' 
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substituting into Eq. (la), and collecting like powers of n yields the follow* 
ing equations for the flapping behavior away from the boundaries. 




[(1 - r2)u^l* + 2 u) 2 u^ ■ 0 s Xq 

[(1 - r2)u|]‘ ♦ 2(0^ > -2«2 = Xj 

[(1 - r2)u[]' + 2«J u^ » - 2 jCu? u^_^ = x^ i » 2,3,. 


( 6 ) 


Eq. (6) is simply the well-known Legendre equation with a nonzero right-hand 
side. Its solution can be expressed as integrals of the Legendre functions, 
P 5 


V • 


U. a C P ♦ P 

1 1 V V 


So ' 




2ui^ ■ v(v ♦ 1) V > 0 

o * 


The constants Cj, C 2 , and v of Eq. (7) must be determined from the bound- 
ary conditions of the problem. The condition Ui finite on the interval 
0 <r < 1 requires: (1) that v be an integer tin which case the reduce 

to Legendre polynomials), and (2) that 


/ I 

p X. d; 

V ^1 

o 

(so that the indefinite integral remains bounded). The boundary condition 
® requires that v assume only odd values, because the even Legendre 
polynomials are nonzero at r « 0. Finally, the constant must be chosen 
so that the order 1 solution is not repeated in the higher order solutions, 
thus avoiding secular terms. Therefore, Cj ■ 1 for i ■ 0 and C » 0 for 
i 0. The complete solutions for all orders of the outer expansion are 
correspondingly 
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il-.ll 


i i 1 I t I I I I i l-i-1 1 I I I - !. i i).. I 1 - 1- 1- i -i;. 


■ n(2n - 1) 


u ■ P 
0 2n-l 


b)| as yet undetermined 


■V. / 


X'. 






dr 


( 8 ) 


n * • • i * • • 

2 

where n is the number of the desired mode, and the u. are the yet undetermined 
frequency expansion coefficients. ^ 

Specific expressions . The first order approximation for the natural fre- 
quencies and modes of a low stiffness beam are given from Eq. (8) as 


mode 

1 

u 


0 

0 

1 

1 

r 

2 

6 

y(Sr^ - 3r) 

n 

n(2n - 1) i 

P , 
2 n-l 


The higher order corrections to u(ui, U 2 , . . .) come from repeated applica- 
tions of Eq. (8) using the properties of Legendre polynomials given in 
Table II. For example, for the first mode. 


^1 * T “i ^ ♦ *“)] 

^2 “ j (- 1 ♦ rtn(l ♦ r)J ♦ j wj dc j tn(l + r) 

j rtn(l + r) ♦ - j iln(2)[2 ♦ r£n(l ♦ r) - r*.n(l - r)]| 


(9a) 


and for the second mode 


Uj • - y ^ r3 - (r - I r3) £n(l ♦ r)j 


(9b) 


(9c) 


In general, the higher order correction terms are combinations of polynomials 
and logarithms with as unknown coefficients; but they do not satisfy the 
boundary conditions. 
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f.' 
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'V- 

If, 

■ 't . 


r* 



1 J ^4 
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In order to define the behavior of (and consequently of u), these 
outer solutions must be matched with inner solutions which themselves 
satisfy the boundary conditions. This matching process will determine the 
behavior of u near r ■ 0 and r ■ 1, and it will also uniquely determine the 
unknown frequency terms It should be noted that the nomenclature "outer" 
and "inner" is mathematical rather than physical. The "outer" solution is 
for the beam region away from the root and tip. The root and tip regions 
each have a distinctive "inner" solution. 

Inner Expansion 

Expansion near root . The yet unfulfilled boundary conditions at the 
beam root are u(0) » 0, u'(0) « 0. From an ordering analysis , the region 
in which these conditions are dominant must be of order To determine 

the behavior of the solution in this innex' region (designated as u to 
distinguish it from tlie outer solution), Eq. (la) must be rewritten in terms 
of an inner variable r ■ 

- 2 + (1 - niP) -2nr~+2nw^u*0 .. (10) 

dP dr dr 

Expanding u in powers of 


u ■ u^ + nUj ♦ . . . 

substituting into Eq. (10), and collecting like powers of n, the equations 
for the uj become 

d**u, d^u 

- 2 =- ♦ i- » 0 


dr 


dP 


d^u^ 

*5^* d?^ 


d*»U2 

d^. - du. , 
p --Hi ^ 2? 1-2 


( 11 ) 


dP 


^ - 2 Z Wj,3 “i. 
j.3 ^ ^ ^ 


dr 


)♦! 


i ■ 3,4,. . . 

The first two expansion variables u, and u, can be quickly found from 
Eq. (11). 2 
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J, . . b,r ♦ c, 


»2 • *j ♦ bjt * C; 


The boundary conditions u(0) • u* (0) • 0 and u remaining finite in the far 
field immediately determine all but one of the constants for each variable. 


“i * N[- 

“2 ' 


(12) 


Thus, the unknown constants to be determined are the b^ in the inner expan* 
Sion and the w| in the outer expansion. These must be chosen so that the 
slopes and displacements of the inner and outer solutions match. 

Expansion near tip . In addition to the geometric boundary conditions 
at the beam root, there are two natural b<Htndaiy conditions at the beam tip, 
u'*(l) ■ u"'(l) > 0, which must be satisfied. An ordering analysis reveals 
that the region in which these conditions are dominant is of order . 
Rewriting Eq. (la) in terms of an inner variable t » (1 - r)/nV^ (and 
designating the solution as 0) yields 


— - f(l - i ~ - (1 - S - n‘/Va - 0 (13) 

df** ^ df2 " 

The variable (I must be expanded in powers of because of the combination 
of with the expansion. 

a ■ 0^ ♦ Oj ♦ Oj ♦ . . . 

Substituting the above into Eq. (13), collecting like powers of n, and taking 
the first integral of the equations yields 
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dQ. 

df5 At 


^ *1 - “l.J * “'f “I /' “= 

j»0 


The boundary condition 


1 d‘’a 


(t - 0) « 0 


implies that the constants of integration a^ must vanish for all orders of n. 

The general homogeneous solution to Eq. (14) can now be expressed in 
terms of Airy's functions. 


■ "i [l -j' >“] ^ J* * “li 


The constants Cj must be zero so that Oj remains finite in the outer field, 
and the constants b^ must be chosen so that the boundary condition 


(f » 0 ) « 0 


is fulfilled. The only available coefficients for matching with the outer 
solutions are simply the additive constants dj. Therefore, only the tip 
displacement for u can be matched, and derivatives of u m <st match automati- 
cally from application of Eq. (14). 


Matching near root , 
formally expressed as 


Matched Expansion 

The matching condition near the beam root can be 


““ u(r) . u(r) 

r -► 0 
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In orter words, if the outer_expansion u is written terms of the inner 
variable r, the expansions u and u mustjnatch term for term to each order of 
T). Now u(^can be expressed for small r as 

u(r) * ♦ u^(0)r] + n[u 2 ( 0 ) + u{(0)r] + . . . 

Comparing this with Eq. (12), as r becomes large, provides the matching 
relations 




uj.,(0) 


(matching elopes) 


u^(0) = - /f bj 


(matching displacements) 


Using the relations in Table II and Appendix A for u(0) and u’(0), the 
unknown constants b^, and b^ can be found for the nth mode. 


b 


1 


u;(0) 


(-1)"'^ (2n)i 
(n-l)J (n)! 2^”'^ 


IlSdi 
/T L 


(2n)! 1 1 

Rni - ijT 


b, = 0 

The inner solution, written in terms of the outer variable r, then becomes 


u = 


—(zll 

(n - 1) 


":!.C2n21 [ 

!(n)! 2^"“^ L 




5"ni/2 
^ - 1 ) 


+ ... 


Since u has been determined for the first mode, for n » 1 may now also 
be found. 

U2(0) “ “ j «2 j “l [j - j in(2)J s b^ ■ 0 


0)2 = 7 - 6tn (2) » 2.84 , (n » 1) 


10 




Matching the root solutions has therefore given the ““>de and frequency 

to order n and higher modes and frequencies to order 

Matching near tip. The matching conditions near the root have specified 
0)2 and 0 ) 2 . ^e next step is to match the tip expansion to the outer expan- 
sion and^thus determine the effect of tip boundary conditions on the frequency. 
Writing u in terms of t yields for small t 

u - Uq( 1) ♦ + n2/6[-u;(l)f] ♦ n3/6[u^(l)] + (Dr^] 

t nSZ6[- u}(i)t] + |u*"(1)t3 ^ + . . . 


Repeated applications of Eq. (14) yield 


^0 ' *^0 “ ^ 


0 ad 


Qn ® d 

2 2 


0)^ f « - wj f 
0 o 




\ “ I "o^“2 ■ 


[t* AjW [3 


Q 5 ’ - f[“l + “o u^(l)l. d, 


w^(u)^ - l)(u^ - 
0 0 ' ' r\ 


2 - 3) r.3 ft 

^ f Gi dr 

/ •'o J 


‘‘e - 


+ integral terms 

These solutions automatically match the higher powers of n, for example 

u;(l) » u»2 

u'*(i) -io,2(u)J-l) 

u'''(l) • |u.2(u)J - D(u)J - 3) 

Uj(l) - 0)2 + o)2Uj (1) 
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Lij f n n f 1 1 1 r-n'' m 



so that the only consttmts required for amtchlntt »ve the previously vletevmlneU 


The solutions provide some useful lnforB\»tlon voneernlnii the tip region. 
First, while fulfilling u*lO) ■ i' introduced expotientlals, fxilfilllng 
u»'(l) » 0 and u***(l) » 0 has introduced Ai and di into the solution. 

Second, although Ai does vlecay e,\ponentiallv Ifor large fl dl only decays as 
l/t (for large t), so^that the tip expansion will affect the root slope to 
an algebraic order for n ^ 2 and n*'^'*’for n » 1 (.because u'* - l » o for 
n » 1). The effect of the tip conditions v>n fi'equency is therefore delayed 
into higher oi\lers of n and can be neglected to the order considered here. 

Combin ed express ions . Now that the inner expansions and outer expan 
s i onsTiav e "Been mat ctte4 , they can be combined to form uniform expressions, 
by adding u, u, and d and subtracting like terms, the movlal deflections (and 
frequencies) may be expressed. 

r. ■ I 


I * ♦ I* - 0 tn(-Ht) ♦ 


n > 2 


. 1) . „1/-T— li'il- - 

■ - 1(|0U1> ■ Ul 


->.-T 

>,'n-lj 


♦ iH'i) 


n » I 


u “ r ♦ »'T * oxp(- r/*T 

* ' 1*^/* ' ♦ r - 1'^ - ♦ ri 

^ 0 

'] 


-I’ - <nl^))r <nO ♦ r)| ♦ Oitr*^-’) 


(ISl 


riviUl * il 


.. l..^ t*-'* ' 1/ ' r T '' ' , S ,s 

u - vlJ“ - •'«) ♦ '!*'■ j vj • ‘ ♦ »■ - -.n* ■ ^ «■ ) 

- Vj- expi- r/*'Jn*'^')j ‘ (rl - i:urlA|(rlj 


‘ 0 Ml 


A comparison of the asymptotic fi-equency expression in Iq. (ISl with 
dalerkin's method results Is given In Fig. 1‘his comparison reveals 


that below » 0.15 the n ■ 10 Galerkin results break down. The figure also 
reveals that the 0(n) expansion is within 10% of the actual value even up to 
“NR • 0«45. It may be concluded that for the first mode the expansions show 
good convergence. In Fig. 3, a comparison of Galerkin and expansion mode 
shapes is given. For wjsjg * 0, the Galerkin results naturally fail to predict 
the actual straight line mode shape. For wfjR » 0.35, in which the Galerkin 
results can be assumed to be nearly exact, the asymptotic expansion gives 
an excellent approximation even to order . Therefore, the first mode 
shape also shows good convergence. 

Two valuable pieces of information result from this asymptotic expan- 
sion. First of all, Eq. (15) gives the behavior of and u for small n. 

This is necessary information for obtaining a simple frequency expression 
valid at all n. Second, the form of the inner expansion can be used to 
estimate the required number of orthogonal functions needed for convergence 
of the eigenvalue determinant. Since the expansion indicates an 




it follows that 


Sm > ^ 


for convergence. Since = tt(N - 1/2), a simple expression for the required 
number of modal functions can be obtained. 


N > i = j * — 

2 L J 2 


Figure 1 confirms the validity of this formula. 


EMPIRICAL COMPOSITE EXPRESSIONS FOR BENDING 

Two steps now remain in determining a simple frequency expression for 
flapping. First, the behavior of for large n must be found. The method 
used here is to expand the eigenvalue determinant in powers of n. Second, 
a search must be made for a simple empirical function which can model both 
the low n and high n characteristics. 

Large Stiffness Expansion 

Eigenvalue determinant . The determinant of Eq. (3) is characterired 
by the tact that, for large n, it is dominated by its diagonal elements 

D.^(i s j) » D^^(i ^ j) 

In this case, the eigenvalues of D (A) may be written as 


n nn n 


6A « A 

n n 
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i; I i\i t r t i' 


.t. j.vH- f -iViI’ 


The eigenvalue determinant, therefore, can be rewritten in terms of 5A 


where 


[hi] - y'lj = Kj] ‘ “ni'i • « 


hi ‘hi 


D. . * D. . - D i * j 
13 xj nn •’ 


'y • ») • 


Expanding the determinant, k( , ing only the highest power of D. .(i = j) in 
eiv h ten" , gives ^ 


(5A + (6A D.. ♦ 

n^ n-' 11 

i/n 


* °ii 




But 6A is of the order of the off diagonal terms, so that to first order 
only the last two t rm£ in E£. (16) need be retained. Discarding the higher 
order terras and using 


ni' 11 

1*1 

ij*n 


, _ n V* D'-ni 

” ' ' i?i “ii ■ ”, 

i/n 


The eigenvalue, therefore, is given as a summation of combinations of the 
elements of D. 

Expansion tor flapping . When the are taken to be nonrotating modes, 
the frequency expansion for flapping becomes 


(.32 = A 


. ne; - 2: 

it 
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yi iigi 


M 


I 


There are two types of truncations associated with Eq. (18). The first 
type of truncation, from Eq. (16), is of order less than (1/n)^. The 
second type is the truncation of the summation in Eq. (18). For example, 
the frequencies for n - 1,2 from Eq. (18) are 


n a 1 


= 1.19 ♦ 12.36n 


.0890 


.0376 


.0085 


1 + 89. 3n 1 + 228. n 1 + 419. n 


n a 2 


a 6.48 + 485. n + 


.0890 


.0025 


.2860 


1 ♦ 89. 3n 1 + 292. n 1 ♦ 476.n 


Truncating the series is similar to truncating the number of orthogonal func- 
tions. The truncation error will clearly be much smaller than 1/n, as it is 
bounded at a fixed percentage as n goes to zero. The formulas in Eq. (19) 
are compared, in Fig. (4) and Fig. (5), with the exact eigenvalues of D for 
several truncation values N. The results indicate that Eq. (19) is a good 
approximation for the eigenvalues of D even down to very small n. 

Similarly, approximate eigenvectors can be determined by 


[[D] - A^[I]] {a^} = [D){a^} = {0} 


from which it immediately follows that for flapping 


A.. - A +n(6'? -6 ) 
11 nn ^1 n-* 


i I* n 


Eq. (20) also agrees rather well with the Galerkin results. For example, 
comparing the expansion modes from Eq. (20) 


.130/(1 + 89. 3n) 


(-.130/(1 + 89. 3n) 


.048/(1 + 228. n) 


-.015/(1 + 292. n) 


.002/(1 + 419. n) 


.099/(1 + 476. n) 


with the Galerkin modes for n = 0 
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-J 



il 






it is clear that the modal expansion converges very rapidly. Therefore, 
Eq. (18) and Eq. (20) are seen to be excellent approximations to and 


a^j for moderate to large n. 


Approximate Frequencies 

Choice of function. Now that the behavior of is known for small and 


large values of n> it remains to find a function which will be valid at all 
values of n. When the order of the infinite series in Eq. (18) is estimated 
by using an integral rather than a summation, products of and tan~^(ri^'^) 
are introduced. These possess potential for use in a composite function for 




Of particular interest is the combination 


s a bn + cn^^^ tan~^ (d/n^/^) 


( 22 ) 


This function has the advantage of possessing an expansion in n^/'^ for small 
n and an expansion in n for large n. The next step, therefore, is to deter- 
mine if Eq. (22) can be made to approximate the behavior of 


Choice of constants. For small n. Eq. (22) can be expanded as 


„2 B a ♦ ♦ (b - c/d)n ♦ . . . 


(23a) 


To fit the first two terms to the small ti asymptotic expansion, Eq. (15), 
requires 

— - r 

a = n(2n - 1), c « ^ ‘ 


r {2n)! 1 1 

ir [(n)j(n - l)j 


For large n> Eq. (22) can be expanded as 


* bn ♦ (a + cd) - 


cd^ 


n 


(23b) 


To fit the first two terms to the large n eigenvalue expansion, b and d 
must satisfy 


S' 


nn 


n(2n - 1) 


In terms of » 3^ n^^^ Eq. (22) then becomes 


J 
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n(2n - 1) 
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nn 


n(2n 


- 1 ) 




(24) 


c/ei 


Thus, Eq. (24) represents an approximation valid for all frequencies. The 
terms in this approximation are simply ujjr (the nonrotating frequency), 
n(2n - 1) (the small n centrifugal stiffness), (the large n centrifugal 
stiffness, and f (the order frequency correction). The equation applies 
equally to lead-lag with replaced by 


Evaluation of results . Eq. (24) is a simple empirical expression rep- 
resenting 0)2 and is compared with Galerkin's results (for n = 1, 2) in 
Fig. (6) and Fig. (7). It is clearly seen that Eq. (24) is a good fit to 
0 )^ at all n values. This closeness of fit may also be viewed quantitatively 
by comparing the third terms in Eq. (23a) and Eq. (23b) with the corrc-spond- 
ing terms in the asymptotic and eigenvalue expansions. For example, the 
order n term, for small values of n and n = 1, is found to be 


b - c/d * 12.36 - 9.44 « 2.92 (composite) 

7-6 in(2) = 2.84 (asymptotic) 

Similarly, the order 1/n term for large values of n is found to be 

- cd^/3 = -.0013 (composite) 

.089 .0376 .0085 


= -.0012 


(eigenvalue) 


89.3 228. 419. 

Eq. (24) is therefore an excellent fit to the exact frequency behavior. 

Approximate Modes 


Observations . The final step in solving the flap equation is a deter- 
mination of approximate expressions for the mode shapes u(r). In Fig. 3, 
the low stiffness first mode (to order ^) and the high stiffness first 
mode (to order l/n) are seen to be virtually identical. Similarly, Fig. 8 
provides a comparison of the second mode. Except near the blade root, the 
high stiffness expansion (to order 1/n) is very good even for n = 0. 
Furthermore, the mode experiences only minor changes as n is increased. 
Because of the n range where the two modes are in close proximity, it is 
possible to form approximate modes as linear combinations of the two expan- 
sions. 

Normalization . Before combining the small and large n expressions, it 
is necessary to normalize the modes such that 

r u^ dr ® 1 

*0 
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In Eq. (20) and Eq. (21), the high stiffness modes are seen to be completely 
dominated by the nth nonrotating mode (which is itself normalized) so that 
no normalization procedure is necessary. The small stiffness modes, however, 
are not normalized so that some additional work is required. 

For all but the first mode, only the first order mode shape P 2 n-i 
be used. Table II provides the relation for normalizing these modes. 

n > 2 


u - y4n - 1 (-1)"*^ (2Sa) 

(The negative sign is added to give the small stiffness modes the same sign 
as the nonrotating modes.) For the first mode, the expansion will be taken 
to order The normalized version of this mode, neglec* higher order 

terms, can be calculated to be 


n e 1 


u ■ vT |r ♦ 1 ^ j " 2 in (2)^ r ♦ rUn (1 + r) 

+ exp(- t//2t\ 


(2Sb) 


Eq. (20), Eq. (2Sa) , and Eq. (2Sb) are the normalized modes which will be 
combined to form the composite modes. 


Modal composition . When two normalized modes are nearly the same, having 
two si mi 1 ar frequencies , the formation of a third normalized mode having some 
intermediate frequency can be accomplished by a simple linear composition. 
Designating w_(Ug) ao the smaller frequency (with its mode) and w, (u,) as the 
larger frequency^ (and its mode), the desired mode u having frequency u is 
reasonably chosen as 


u « u 


S 




(26) 


As u approaches either frequency, the mode corresponding to that frequency 
dominates. 


In the context of this paper, Ug is taken to be the mode of Eq. 
and u)g is taken to be the asymptotic^expansion frequency. 


( 25 ), 


n ■ 1, 
n > 2, 


“S 

.,2 


J + 1 . ,,1/2 + [7 - 6 «n (2)]n 

/2 

„(2„ . i) . ivi 
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Similarly 0 )^ and ul are taken from Eq. (18) and Eq, (20). The choice of 
as given by Eq. (24) insures that w| < making the linear composition 

of Eq. (26) a smooth transition from Legendre polynomial (n ■ 0) to nonrota- 
ting mode (n ■ »). It should be pointed out that by replacing with 

.Ar. ♦ 1* these modes also apply to the lead-lag case (see Eq. (1)). 

LcAD-LAu 


ANALYTIC SOLUTION FOR TORSION 

Now that the flapping and lead-lag solutions are in hand, the solution 
of the torsion equation remains. The torsion equation from Eq. (Ic) is 

-ye" - j[(l - r2)0»]' - (u)2 - 1)6 « 0 

The first term, the torsional rigidity term, is only of second order dis- 
tinguishing this equation from the fourth order bending equations. Tlje 
second term is the rigidity from centrifugal tension, and the -1 is the 
so-called "tennis racket" term resulting from gyroscopic forces in the 
rotating reference frame. Although the torsional frequency (and the stiff- 
ness y) are often large enough that effects of rotation can be neglected, 
the lower limit of practical frequencies does include significant departures 
from the nonrotating case. The task at hand is to develop simple frequency 
and mode shape approximations which will: 1.) be valid for all values of 

torsional stiffness and 2.) thereby determine the range of y and u for 
which rotational effects may be neglected. 

Closed Form Modes 

Limiting cases. Before proceeding to the direct solution of the torsion 
equation, it is interesting to study its behavior in the limits of y = 0 and 
y n ®. When y is identically equal to zero, the solutions to the torsion 
equation are Legendre polynomials. 

* 1 + n(2n - 1) 

6 “ ^ 1.2,3,. . . 


The lower limit, therefore, for the torsional frequency of uniform rotating 
beams is In the limit as y becomes very large, the torsion solutions 

become the nonrotating modes and frequencies given in Table II. 


? _ (2n - 1 ) 2^2 

' • 4 > 





These nonrotating modes, unlike the y **0 modes, fulfill the boundary condi- 
tion 6(1) • 0. 

Change of variable . In contrast to the flapping and lead-lag equations, 
the torsion equation aclmits a closed form solution for the mode shapes. By 
a simple change of independent variable 


r ■ xv'l ♦ 2y , X « r//I~+Tv 

the torsion equation and its boundary conditions can be brought into the form 


6 ( 0 ) ■ 0 


de 
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\»^1 + 2 y / 


1) « 0 

(27a) 

0 . 

(27b) 


The torsion equation has, therefore, been transformed into Legendre's equa- 
tion. The beam tip, however, is at x » - rather than x = 1. The cen- 

/l ♦ 2 y 

eral solution to Eq. (27a) is, from Eq. (7), 


6 « aP 

V 


♦ bP 


r 


Pjci - x2) 


dx » aP ♦ bQ 


(28a) 


2u»2 - + 1) + 2 , V > 0 . (28b) 

From the properties of Legendre functions, the boundary condition 
6(0) ■ 0 implies that 


./b . i tan (f ) 

yielding a series representation for 6.® 


[«- 


* ■ 2 ) „3 * (v - 1) (V - 3}(v 4 2)(V ♦ 4) 5 


3! 


S' 




(29) 


The remaining boundary condition, ^ {l/^V ♦ 2y) » 0 is fulfilled only for 

certain discrete values of v. These eigenvalues for v then imply the fre- 
quencies and modes through Eq. (28b) and Eq. (29). 



The mode* of Eo. (29) are valid for all values of v. recovering both 
limiting case* (y ■ 0, •). When y ■ 0, Eq. (27a) becomes singular and the 
boundary condition is relaxed. The series expansion of Eq. (29), however, 
diverges unless v takes on one of the eigenvalues 

V ■ 2n - 1 n ■ 1,2,3 


This results In a truncation of Eq. (29) yielding the Legendre polynomials. 
For large values of y, Eq, (29) reduces to 

0 ■ ^ sin^vr/*^ Y^^^J 

The boundary condition 0*(1) ■ 0 then implies 



(2n - l)n 

■7 



(2n - l)^Tt" 

4 


Y 


recovering the nonrotating solution. 

Approximate Frequencies 

Small Y expansion . Eq. (28) and Eq. (29) now describe the torsion mode 
in closed form, out the solution of a transcendental equation 



is required to determine the natural frequencies and undefined constants of 
the modes. In order to obtain a simple frequency expression, this transcen- 
dental equation needs to be approximated in closed form. To begin, the 
behavior of Eq. (27b) is studied for small values of y. 

In the limit of small y, the boundary condition of Eq. (27b) must be 
satisfied for x very nearly 1. 


^TiP * ' ■ 1 - Y ♦ - • • • (^Oa) 

TIP 

Correspondingly, v is expected to be very near to one of its odd integer 
values 

V ■ (2n -l)+e , e<<l . (30b) 

The tip boundary condition, therefore, may be written to first order in , 
using Eq. (28) and Eq. (30), as 
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Using the properties of Legendre functions, this boundary condition becomes 


r , T 

. 2_.l 

■ - tan » ■: p-J 

L jY-‘f(2n-l)inY + 0(l) 


■ ^ tan“^ [y’f{n(2n - 1) [1 n(2n - l)YinY ♦ 0(y)] * 0(c)}] 


■ 2n(2n - 1) (y ♦ n(2n - Dy^^hy] + O(y^) 

Thus, the behavior of v and for small y is found as an expansion in powers 
of Y and YJinY. 


u)2 . 1 ♦ v X v , » , lj , j ^ jj(2n - 1) ♦ n(2n - 1) (4n - 1 )y 


♦ n^(2n • l)^(4n - l)Y^tnY ♦ O(y^) (31) 

^irical frequency expression . The behavior of for large y can be 
found by applying tq. (it) to the torsional determinant, and using Table I 
for determining 

■ I* (2n - 1)2(1 ♦ 3y) 

(32) 

- 192 (2n - 1)** £ (2i - 1)** 

ff2(i ♦ 3 y) 1“1 ((2i - 1)2 . (2n - 1 ) 2 ] 5 

il'n 

As before the summation converges rapidly giving a very accurate approxima- 
tion for u)2. The form of the asymptotic expansion and of the eigenvalue 
expansion provide motivation for attempting to fit the frequency behavior 
with a function of the form 


u)2 ■ a ♦ by ♦ cy2 In 


- e 


d ♦ Y“ 


1 ♦ 3y 


For small y« this formula reduces to 

u>2 ■ a ♦ bY ♦ 2cY^inY - (e ♦ c In (d)]y^ ♦ . 


(33) 
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The constants a, b, and c can be inuaediately determined by comparing with 
Eq. (31). 


a ■ 1 ♦ n(2n *• 1) 
b ■ n(2n - l)(4n - 1) 
c ■ n^(2n - l)2(4n - l)/2 


(34c-c) 


For large y, Eq. (33) becomes 

a.b, - ^ ^ 

. (b - |)t (. - cd . f ) - ^ i . 
From Eq. (32), this gives the values of e and d 


% 2 A 

e » 3n(2n - l)(»n - 1) - ~ (2n - 1)^ 

d . 8n(2n •• 1) (2n ♦ 1) - 3 - 2ir^(2n - 1)^ 
6n^(2n- l)2(4n - 1) 

The final formula for torsion is therefore, 

» 1 ♦ n(2n - 1) ♦ n(2n - l)(4n - l)y 

v2 


(34e,d) 


+ jn^(2n - l)^(4n - l)y^in 


3n2 


d ♦ 


(35) 


(2n - 1)2 - 3n(2n - l)(4n - 1) 


1 ♦ 3y 

To complete the solution, the mode shape is then given by Eq. (29) with 


X • 


/8«u 2 , 7 - 1 
2 


✓1 ♦ 2y 

Practical Considerations 

Comparison of methods . A comparison of the various frequency approxi- 
mations reveals three ^acts. First, the Galerkin expression, Eq. (32), with 
only one mode gives frequencies to within S% for y as low as 0.8 (first 
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torsion frequency 2.0) for n ■ 1, 2. Second, the frequency expansion, 

Eq. (31), gives frequencies to within 5% for y as high as 0.8 for n ■ 1,2, 
And third, a frequency expression neglecting the ((1 - r^)9']’ term 

u)2 ■ 1 + ^(2n - 1)^Y 


will introduce errors of the order of 15% at y * .8. The impact of these 
three observations is double-edged. On the one hand, the low y expansion 
(in terms of logarithms) is seen to be an excellent approximation even as 
Y approaches unity. On the other hand, it appears that for y in the prac- 
tical range of interest the Galerkin solution is perfectly adequate so that 
the low Y expansion is actually not required. This should not, however, be 
construed to imply tha' the effects of rotation are small. The importance 

of the tension term j[(l - r2)0']' reveals that rotation is important, but 
may be accounted for^by a single function Galerkin analysis. Similarly, 

Fig. (9) shows that for y > 0.8 the modes have also converged to the nonrota- 
ting modes so that the mode shape, Eq. (29) is also not necessary. 

Special cases . TTiere may be, however, special cases for which Eq. (29) 
and Eq. (35) are necessary. That is, in Eq. (Ic) the area and mass radii of 
gyration (k/^,k){) were assumed equal. In the event that they are unequal, 
the tension term in Eq. (Ic) must be multiplied by k^/k^. Therefore, the 
region of validity of Eq. (32) will become 

> A * 3 k^/kj 

rather than <d > 2. For example, if k^/kM = 2, then the region of validity of 
Eq. (32) is reduced to w > 3.6. It follows that for beams with certain types 
of cross-sections, )iq. (29) and Eq. (35) may be required if 


u)2 < 1 + 3 k^/k^ 


CONCLUSIONS 


By combining a low stiffness expansion and a high stiffness expansion, 
a simple empirical closed form expression has been obtained for the flapping 
and lead-lag bending frequencies of rotating uniform beams. 


« 1 + 0 ), 
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n B mode number 

■ nonrotating frequency (flap or lead-lag) 

■ Galerkin integral (Table I) 


0 B Galerkin constant (Table I) 
n 

This formula is uniformly valid for all values of from zero to infinity. 

A similar expression has been obtained for torsion, although for most prac- 
tical examples the frequencies are adequately estimated by the Galerkin for- 
mula 

‘4oRSION “ 4 * T2 

These frequency expressions, along with the modal equations which have 
been developed, provide a practical alternative to energy methods for natural 
frequency determination. They provide assured convergence throughout the fre 
quency range, require only minimal numerical calculations, and allow rapid 
determination of the mode shapes and beam properties for any desired rotating 
frequencies . 


v^(4n - 


LZ-ii 


(n - 1)! 
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APPENDIX 


BEHAVIOR OF u^ NEAR BEAM ROOT 


Prom Eq. (8), the formula for u^ (of the nth mode) can be written as 


1 1 2n-iy p2 (1 _ J.2) 


From Table II, P can be expressed for small r as 


P- , * a,r + a,r^ 
2n-l 1 3 


where 


_ _ 1 (-l)”'^(2n)I a - 1 (-l)"~^(2n ♦ 2)1 

1 ” 22 n-l (n)I(n - 1)! * 3 " 22 n (n)l(n + l)! 3 

Therefore, near the root section behaves as 

ft —i-- , !i!! 

u, = 2iof (a r ♦ a r^) I ~ d? 

(a^r^ + 2a a r*‘)(l - r^) 

1 id 

» 2uif (a r + a r3) [-^ ♦ 1 - ^ d? 

(4n - l)a^ Lr^ *1 J 




(4n - l)a 


2 '"1* ”3‘ 


t o2H 


w^(n - l)!(n)! 2 


u. (0) » r , uj (0) « 0 

(-l)"(2n)J(4n - 1) 
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TABLE I. 


A. FLAPPING AND LEAD-LAG ujpjp = 


n 


''n 


. 4 

1 

0.73410 

1.87510 

3.51600 

12.3623 

D 

1.01847 

4.69409 

22.0345 

485.518 

B 

0.9992? 

7.85476 

61.69/- 

3806.55 

B 

1.00003 

10.9955 

120.901 

14617.1 

H 

1 1.00000 

~2 ( 2 n- 1 ) 

( 2 p 1 )^ 

I 

rr4 4 


<) = cosh (p r) - cos (p r) -a„|sinh(p r) - sin (p r)| 
' n n n n n n 

(NONROTATING MODE) 


i 

^1 

^i2 

^3 

^i4 

1 

1.19334 

0.685855 

0.792379 

0.546413 


- 

6.47823 

0.169408 

2.91185 


- 

i 

17.8603 

3.27427 

B 

- 

1 

1 

36.0553 


B. TORSION CP 


NR 


(2n-l);r 1/2 


(2n 1)?T , 

(,-„ = 7 2sm| — 2 — *^1 

(NONROTATING MODE) 


A = { 
i| 


|1/4+^(2i-1)^| 
-4(2i 1)^ (2j 1)^ 


|(2i-l)^ - (2j-1) 


2,2 


I - j 


' ■ ) 
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TABLE II. 


n 

^2n-1 

1 

r 

2 

1/2 (5r^ - 3r) 

3 

1/8 (63r® • 70r^ + 15r) 

4 

1/16 (429r^ - 693r® + 315r^ - 36r) 

GENERAL 

1 -V M)" (2 (2n-1kl|! 

CASE 

2^"*’ ^ (k)! r*n-1-k)! (2n-1-2k)» ^ 

k=0 




n 

1 

2 

3 

B 

GENERAL CASE 

'^2n.1 <« 

1 

•3/2 

15/8 

-35/16 


/o'*2n-l‘^‘' 

1/2 

-1/8 



(-1)" '* (2n-2)! 

4" ’ (2n)((n-1)!j^ 


1/3 

1/7 

1/11 

1/15 

1 

4n-1 
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Figure 1. First frequency from Galerkin's method 
















Figure 4, First frequency from large stiffness expansion. 
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Figure 5. Second frequency from large stiffness expansion. 
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Figure 6. First frequency from composite expansion 
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Figure 7, Second frequency from composite expansion 
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Figure 8. Second flap mode shape 
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